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Abstract. We study a Langevin equation describing the Brownian motion of 
an object subjected to a viscous drag, an external constant force, and a solid 
friction force of the Coulomb type. In a previous work [H. Touchette, E. Van der 
Straeten, W. Just, J. Phys. A: Math. Theor. 43, 445002, 2010], we have presented 
the exact solution of the velocity propagator of this equation based on a spectral 
decomposition of the corresponding Fokker-Planck equation. Here, we present 
an alternative, exact solution based on the Laplace transform of this equation, 
which has the advantage of being expressed in closed form. From this solution, 
we also obtain closed-form expressions for the Laplace transform of the velocity 
autocorrelation function and for the power spectrum, i.e., the Fourier transform of 
the autocorrelation function. The behavior of the power spectrum as a function of 
the dry friction force and external forcing shows a clear crossover between stick 
and slip regimes known to occur in the presence of solid friction. 



PACS numbers: 05.40.Jc, 46.55.+d, 05.10.Gg, 02.30.Jr, 02.50.Cw 
1. Introduction 

We continue in this article our study of Brownian motion involving solid (dry or 
Coulomb) friction in addition to viscous friction; see [l]-[3] . As in these works, we 
consider the piecewise linear Langevin equation 

i^-fv-Aa^+F + VT^t), (1) 

where 7 > denotes the viscous coefficient, A > the dry friction coefficient, F 
an external constant forcing, a Gaussian white noise and T the related diffusion 
constant. The term — Aer(v), where cr{v) denotes the sign of v with the convention 
cr(0) = 0, represents the dry friction force. Its physical interpretation follows by 
considering Eq. ([I]) in the deterministic limit T = 0: for |F| < A, the stationary state 
of this equation is the "sticking" state v = 0, whereas for |F| > A, the stationary state 
is a "sliding" state with v ^ 0. Thus, to induce motion (v ^ 0) from rest (v = 0), F 
has to be larger than A, the dry friction contact force||] 

The effect of dry friction on the properties of Brownian motion was studied by 
De Gennes [I], who showed, for the special case 7 = F = 0, that the velocity- velocity 

X Mathematically speaking, Eq. |TJ is incomplete: in order to ensure the existence of well-defined 
global solutions of this equation, we must ask in addition that all external forces vanish in the stick 
state v = when |_F| < A. 
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correlation function (u(t)u(O)} acquires a dependence on the noise power when A > 0. 
In [3|, we extended his study by obtaining eigenfunction expansions of the propagator 
p(v, t\vo,0) for the general case 7 > and F > 0. From these expansions, which 
involve a special function known as the parabolic cylinder function, we were also able 
to compute (v(t)v(0)). Our main finding was that the stick and slip states of the 
deterministic system (L = 0) translate into stick and slip regimes in the noisy system 
(r ^ 0), which are characterized by a strong and weak dependence, respectively, of the 
correlation time of the velocity-velocity correlation function with the external force F. 

Our goal in this paper is to show that our exact results of [3] can be expressed in a 
more convenient way by solving the Fokker-Planck equation associated with Eq. ([I]) in 
Laplace space rather than in direct space. The resulting expressions for the propagator 
and velocity-velocity correlation function are indeed somewhat more compact and more 
elegant than the eigenfunction expansions presented earlier. Our solution in Laplace 
space also enables us to complete the study of Eq. ([T]) by deriving exact, closed-form 
expressions for the power spectrum of this equation, i.e., the Fourier transform of 
the velocity-velocity correlation function. These expressions extend early articles by 
Caughey and Dienes [5] and by Atkinson and Caughey [6j[7], recently brought to our 
attention, which considered the power spectrum of the Langevin equation with pure 
Coulomb friction, i.e., the same case (7 = F = 0) considered by De Gennes. 



2. Propagator 

As in [3], we study the propagator P(x, t\x' 7 0) of Eq. ([!]), expressed in terms of the non- 
dimensional variables ^2j/Tv — > x 7 jt —¥ t. The Fokker-Planck equation governing 
the evolution of this propagator has the form 

dP d d 2 P 

- = + (2) 

where 6 = A^/2/( 7 r) measures the magnitude of the dry friction relative to the viscous 
damping while / = Fj ^2/ (~fT) stands for the external constant force, measured also 
against the viscous damping. 

The time-independent or stationary solution of the Fokker-Planck equation has 
the standard form 

e -*(a0 

Pf(x) = -g- (3) 
in terms of the potential 

*(*, = M+fi!_ /s . w 

Depending on the sign of the parameter 5, the potential refers to a dry friction problem 
(S > 0) or a Kramer-type tunneling problem (5 < 0){f] 

For piecewise-linear Fokker-Planck equations, the propagator can be obtained 
exactly, as was observed in (mm, by considering its Laplace transform: 



P(x,s\x',0)= / e" st P{x,t\x' ,0)dt. (5) 
Jo 

With this transform, Eq. ([2| becomes a second-order ordinary differential equation 

sP(x, s\x', 0) - S(x - x 1 ) = ^-{x + Sa(x) - f)P(x, s\x', 0) + (6) 

ax aar 

§ Although we do not study the case <5 < 0, all results derived here are also valid for this case. 
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Figure 1. (Color online) Propagator P(x,t\x' , 0) in real time obtained by 
numerically inverting the Laplace transform P(x, s\x' , 0). Parameters: 5=1, 
/ = 0.5, x' = 3. The different colored curves moving to the left are obtained for 
increasing times r = 0.25, 0.5, 0.75, 1 and 2. The dash curve corresponds to the 
stationary density Pf(x). 



which can be solved in terms of parabolic cylinder functions. The detail of this solution 
is given in |Appendix A| The final result obtained for positive values of the initial 
condition, i.e., x' > 0, has for expression 



P(x, s\x', 0) = 9< ( S) S, f) e (*'+*-/) 2 /4 D_ s (x' + 6-f) 



/2tt 



for x < 0, 



x e-( x - s - f f/ 4 D- s {-x + 6 + f) 



P(x, s\x', 0) = 9> ( S) S, f) ^ e ^' +s -^' i D. s (x' +6-f) 



/2tt 



for < x < 



and 



xe -( x + s -tf/ 4 D^ s (x + 5-f) 
+ F j*leV +5 -tf^D_ s (x' + 6-f) 

x e- ( * x+& -rfl A D_ s (~x-5 + f) 



P(x, s\x', 0) = 9> (s, S, f) ef+s-W D_ s (x' + S - /) 



/2tt 



x e-( x+5 -V 2 / 4 D_ s {x + S- f) 

+ r M e (^-/f/^_ s (_ x '_ i + / ) 
V 2n 



(7) 



(8) 



x e -b+*-t) / 4 D_ s (x + <5 - /) (9) 

for x > x'. The coefficients g< and g> are defined in Appendix A by Eqs. ( A.21[ ) and 
(A. 23), respectively. The propagator for negative values of the initial condition, i.e., 
x' < 0, is obtained from the equations above simply by replacing x, x' , and / by — x, 
—x', and — /, respectively. 

The set of Eqs. ([7|-(|9| is the main result of this paper. To gain some insight into 
the form of this solution, consider the case without dry friction and external forcing, 
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i.e., S = and / = 0, which corresponds to the Ornstein-Uhlenbeck process. Then 



3< (s, 0, 0) = 1 and g> (s, 0, 0) = 0, according to Eqs. ( |A.21[ )- (|A.23p of |Appendix A| and 
we are led tcjjj] 

P(x six' 0) r(s) I e {x ' )2/i D- s {x>)e-* 2 /*D_ s {-x), x < x> 

which is the Laplace transform of the Ornstein-Uhlenbeck propagator. Comparing 
this particular form of P with the general solution above, we see that the products 
of parabolic cylinder functions occurring in Eqs. ([7|-([9]) are essentially the Laplace 
transform of the Ornstein-Uhlenbeck process, modified to include the dry friction 
and external forces. Following the method of images, the additional term with the 
coefficients g < and g > can then be considered, at a superficial level, as convolution 
integrals in the time domain, with g < and g > playing the role of the source terms. 

From the relatively compact solution of the propagator in Laplace space, we are not 
able to obtain the propagator itself in closed form, except for special cases, such as the 
Ornstein-Uhlenbeck process (5 = f = 0) and the pure dry friction case corresponding to 
7 — F — 0. However, for all cases it is possible to obtain the propagator numerically by 
inverting the Laplace transform. Figure [l] shows the result of this procedure using the 
so-called Talbot method [8-10 ^] This figure reproduces exactly our previous results for 



the propagator based on the spectral decomposition of the Fokker-Planck solution (see 
Fig. 10 of (3]). In general, we have found that the numerical computation of P{x, t\x', 0) 
from the inverse Laplace transform of P(x, s\x' , 0) is stable and can be carried out to 
arbitrary level of accuracy for a large range of physically-relevant parameter values. 
Our Laplace solution can therefore be considered a useful complement to the spectral 
solution presented in |3]. 

3. Power spectrum 

Using the closed-form solution for P(x, s\x' , 0), we now derive an expression for the 
Laplace transform of the auto-correlation function: 

/>oo 

C(s) = / e - st (x(t)x(0)) alt. (11) 
Jo 

In terms of P(x, s\x' , 0), we thus have 

/oo />oo 
dx' I dxxx 1 P(x,s\x',0) p f (x'), (12) 
-OO J —OO 

where pf(x) is the stationary density of the Fokker-Planck equation, given by 

l f e -(-+S-f) 2 /2 e (S-ff/2 if X>Q 
Pf[X) ~Z \ e -(-*-/) 2 /2) e (S+ff/2 if x < {U) 

with normalization 

D-i{6 + f) D-i(6-f) 

D (S + f) + D Q (S-f) ■ [ > 

|| The steps leading to this result yields some interesting integral identities for parabolic cylinder 
functions, which, to the best of our knowledge, cannot be found in the literature. 
K A Mathcmatica implementation of this method is available at 
http : //library . wolf ram. com/inf ocenter/MathSource/5026/ 
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This expression of Pf{x) is invariant with respect to the inversions x — > —x and / — > — /, 
as is the propagator. As a result, we can rewrite C(s) as 

/■00 /"OO 

C{s)= dx' dx(xx' P(x,s\x',0)p f (x') + 'f -> -/'), (15) 

JO J — oo 

where the symbol '/ — > — /' indicates the contribution to the kernel obtained by 
replacing / by -/. 



We proceed to evaluate the integrals in Eq. (15). The first one in x can be 



performed using a known relation for the derivative of parabolic cylinder functions, 



which can be obtained from the relation shown in Eq. (A. 17) of Appendix A The 
result is 



dx' 



dx x x P(x, s\x' , O)pf(x') = / dx' 

Jo 



(x'f x'(x' + S-f) 
s s(s + 1) 



D- s - 2 {6 - f)T(s) 
Zx/2tt 



e -(S+f) 2 /4 
e -(<5-/) 2 /4 



g< (s,6,f)D_ s _ 2 (6 + f) 



+ g> (s,6, f)D_ s _ 2 (8 - /) + D_ s _ 2 (-<5 + /) 



(16) 



The second term on the right-hand side of this expression can be simplified using the 
definitions ( |A.21 ) and (A.23) as well as the identities (A. 17 1 and (A. 19): 

/ e -(<5+/) 2 /4 

D- S - 2 (S - /) - (s _ f)2/4 9<(s, S, f)D-.- 2 (8 + f) 



+ 9> (s, 6, f)D_ s _ 2 (S -f) + D_ s _ 2 (-6 + f) 



t(* + /) 



s(s + l)r(s) D_ s (6 - fiD.^S + /) + D_ S (S + /)£)_ s _ x (5 - /) ' 



(17) 



By combining this result in Eq.(15), and by performing the integral over x', we then 
obtain 

C(s) = 1 (x 2 ) f - - 7 -i^ ((x 2 ) f + 8{\x\) f - f(x) 



25 



s(s + l) VN ' 
1 £>_ s _ 1 (<5 + f)D 



H\ x \)f - j \x/f) 

2 (8-f)+D. 



- f)D-.- 2 {8 + f) 



1 s(s + l)Z D- a {8 + f)D- B - 1 {8-f) + D- a {8-f)D- a .. l {5 + f) ' (18) 
where (•)/ denotes the expected value with respect to the stationary distribution pf. 
Surprisingly, all of the stationary expected values appearing above can be conveniently 
written in terms of parabolic cylinder functions: 

MS + f)D_ 3 (8 -f)+ D (8 - f)D_ 3 (S + f) 



(N) 



D (S + f)D_ x (8 -f)+ D Q (5 - + /) 

= Do(5 + f)D- 2 (S - f) + D (6 - f)D- 2 (5 + /) 
; D (5 + - /) + D (.5 - f)D-x{S + f) 

D (8 + f)D- 2 (S - f) - Dq(5 - f)D- 2 (8 + /) 
D (6 + f)D_ 1 (8 -f) + D (8 - f)D^[S + /) ' 



(19) 
(20) 
(21) 



The expression shown in Eq. ( 18 1 can be re-written in other ways to make some of 



its properties more explicit. Using the identity (A. 17), it is easy to check, in particular, 

(22) 



that the stationary expectation values obey the relation 



(x i )f + S{\x\) f -f{x) f = l, 
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so that 

p/ \ = \ ( i 2\ _ ! , J_ s _!(^ + f)D- a - 2 (S - f) + £>- a -i(<* - f)D-s-2(S + /) 
lSj a ^ ;/ £ + /)I>- a -i (<* - /) + - f)D- s -!(S + f) 

25 D- B -x(8 + f)D_ s _ 2 (S - /) + D_ s _ x (6 - f)D_ s _ 2 {8 + /) 



IV 2 + /)£>- s -i (* -/)+-»-. (* ~f)D- s -i (S + f) )■ (23) 

This form of the Laplace transform of the correlation function, or resolvent, is 
quite useful to uncover its analytical structure. Because of Eq. ( 14 ) and the identity 
(A. 17) (for v — 0), the coefficient of l/(s + 1) vanishes at s = — 1 for <5 7^ 0. Thus, 
the apparent singularity at s = — 1 is removable. For a similar reason, the coefficient 
of 1/s at s = 0, which corresponds to the residue at s = 0, is given by (x)f. Thus, 
the long-time limit of the correlation function is given, as expected, by the square of 
the mean velocity. All the other poles of the resolvent C(s) are determined by the 
denominator of the expression above, i.e., by the zeros of 

D_ S (S + f)D- a -!(6 - /) + D_ a (S - + /), (24) 

which is precisely the characteristic equation derived from the Fokker-Planck operator 

i- - 

Another form of C(s) can be obtained by rearranging Eq. (23) using the 
aforementioned identities to obtain 



C(s) 



(x 2 ) f - 1 1 2S/Z 



s s+1 s(s + l) 2 

45 2 /Z ( D_ s (8 + f) D-.(6-f) 



s(s + l) 2 {D—itf + f) D_,_ x (*-/) 



(25) 



This expression is better suited for numerical calculations than either Eq. (21) or (23). 
The limit S = f = is also clearer at the level of this expression. Noting from Eq. (22) 
that (x 2 ) f — 1 in this limit, we recover C(s) = l/(s + 1), which characterizes the 
simple exponential decay of the correlation function of the Ornstein-Uhlenbeck process. 
Finally, we obtain from Eq. ( 25 1 a simple analytic expression for the power spectrum 
p(oj) = ReC(s = iuj), i.e., for the Fourier transform of the auto-correlation function, 
namely, 

1 45 /Z 



l+LU 2 (1 

4<5 2 /Z T / 1 / D_fa,(<y + /) D- iu {5-f) 



p(u) = 



uj Im ^(l + iu)2 U— 1(<5 + /) + D^ 1 (5-f)J J' (26) 
This result is used in the next section to discuss the stick-slip transition occurring at 



4. Stick-slip transition 

We have discussed in detail the behavior of the correlation function (x(t)x(O)) as a 
function of 5 and / in [3j and, more precisely, how the stick-slip transition that appears 
in the deterministic (r = 0) equation when / = 6 is modified in the presence of noise 
(r > 0) to a smooth crossover between a stick and a slip regimes, characterized by 
different exponential decay of (x(t)x(0)). Figure [2] shows how this crossover shows 
up at the level of the power spectrum. We see that p(u) is rather flat in the stick 
regime (/ < <5), and that it starts to develop a sharp zero- frequency peak at the 
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Figure 2. (Color online) 3D plot of the power spectrum p(oj) normalized by the 
variance Ax 2 obtained for <5 = 12 as a function of frequency u> and the external 
forcing /. The stick-slip transition line / = S is shown as the red line. 
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5 10 15 

Li 

Figure 3. (Color online) Solid lines: Log-linear plot of the power spectrum p(u>) 
as function of the frequency u for 8 = 12 and different values of the external forcing 
/ (shown in th e pl ot). Colored dashed lines: Approximation of the power spectrum 
given by Eq. \29) . Top dashed line in black: Power spectrum of the Ornstein- 
Uhlenbeck process without dry friction. Inset: log-log plot of p(co) showing the 
uj~ 2 tail behavior. 




stick-slip transition / = 5. This peak is the translation in frequency of the increase 
of the correlation time associated with (x(i)x(0)) as we go from the slip to the stick 
regimes j3j. 

To gain further insight into the behavior of p(ui) , we plot this function on a log- 
linear scale in Fig. [3] together with an asymptotic expansion of this function obtained 
in the limit of large dry friction 5 and large forcing /. Mathematically this expansion 
is equivalent to the smal l no ise limit of the Langevin equation, and is obtained by 



using the representation (23) for C(s) to express the power spectrum as 



p(oj) = — Im 



uj \ 1 + iui 

' 26 F 1+w (5 + f) + F 1+l0J (S - f) 

*F (5 + f)+ F Q (S - f) 26 + (1 + ilo){F 1+ilc {6 + /) + F 1+iui (6 - /)) 



(27) 
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Figure 4. (Color online) Left: Simulated paths of Brownian motion with dry 
friction (for the rescaled variable x) for <5 = 12 and / = 10 (blue, stick regime) 
and / = 14 (purple, slip regime). The numerical integration was done with the 
Euler-Maruyama method with At = 0.01. The steady-state velocity for the slip 
regime is x* = f — 5 = 2. Right: Corresponding power spectra (colored lines) 
compared with the theory (black line). The numerical spectra were obtained by 
averaging 50 times series over the time interval [0, 100]. 



where 



With some results of asymptotic analysis, presented in | Appendix B we can then obtain 



the following approximation, which is valid for large values of the parameters S and /: 

f il m ( s'-f* 7^+1+ ^+7^+1+^ _ 0<f<S 
PH« -{(^r 2 ^7E^ + ^7E^ +s u -/ <d (29) 

{ V(l +oj 2 ) 0<S<f; 

(see 11 for related results derived via continued fraction expansions). Figure [3] shows 
that this asymptotic formula is relatively accurate, even for rather small parameter 
values. Larger deviations are visible close to the transition point S — f. The inset of 
Fig. [3] also shows that the tail of p(co) decays at large frequencies as w~ 2 , which is 
the sign that temporal correlations decay exponentially, as in the Ornstein-Uhlenbeck 
process. 

To give an idea of how these results might compare in practice with experimental 
results, we show in Fig. [4] simulated paths of the Brownian motion equation with 
dry friction for the rescaled variable x(t) together with their corresponding power 
spectrum. Two paths are shown on the left-hand side of Fig. [4] one in the stick regime 
(S = 12 and / = 10) and one in the slip regime (8 = 12 and / = 14). The power 
spectrum characterizing each of these regimes, shown on the right-hand side of Fig. |3J 
is computed numerically by averaging the spectra of many random paths (here 50) 
over a relatively long time (here T = 100). The numerical results compare well with 
the theory, as can be seen. Instead of averaging different spectra, one can also use, as 
is well known, a frequency window larger than the frequency spacing Auj = 2tt/T to 
obtain a relatively smooth spectrum. 

To close this section, let us now briefly discuss the behavior of two quantities 
derived from the power spectrum. The first is the stationary variance Ax 2 , which is 
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Figure 5. (Color online) Left: Covariance Ax 2 given by Eqs. \l9\ and ( | 2 1 [ ) 
as a function of the dry friction and driving forces. Right: Diffusion constant 
D = p(uj — > 0) as a function of 8 for various values of /. The solid line is the exact 
expression shown in Eq. ( |27| , whereas the dashed line is the diffusion constant 
obtained from the approximation shown in Eq. ||29|l. 



proportional to the total spectral weight: 

, 1 



If 



p(uj)dw. 



(30) 



As shown in Fig. [5] there is a sharp increase of Ax 2 at the stick-slip transition / = 5, 
separating a low variance (stick) regime from a high variance (slip) regime. The second 
quantity is the (non-dimensional) diffusion constant: 



D = 



(x(t)x(O)) dt = lim p(ui). 

uj— >0 



(31) 



This quantity is of particular interest, since it can be measured in experiments. Its 
expression is obtained from Eqs. (27) or (29) and is plotted in Fig. [5} We see that in 
the slip regime, the diffusion constant does not depend much on the external force, as 
is the case for the Ornstein-Uhlenbeck process. This is consistent with the observation 
that the slip regime is essentially a regime of normal Brownian motion in which dry 
friction force plays little role; see [3] for more details. In the stick regime, on the 
other hand, the diffusion constant is relatively small, and sharply increases when the 
stick-slip transition is approached. This behavior is also seen if we derive D from the 
asymptotic expression ( 29 ) . This is illustrated with the dashed lines in Fig. [5] 



5. Conclusion 



We have presented the exact solution of the propagator of a Langevin equation modeling 
Brownian motion in the presence of solid friction, viscous damping, and an external 
constant force. The main feature of this equation is that it shows a stick-slip transition 
often encountered in real systems involving solid friction and external forcing. The 
solution of this equation follows by considering the Laplace transform of its associated 
Fokker-Planck equation, and serves as a complement to a previous exact solution, 
derived in [3] using the spectral decomposition of the Fokker-Planck operator. It also 
extends a similar Laplace solution, previously derived in [5]-[7] for the special case where 
only solid friction is present. 

A clear advantage of the Laplace solution over the spectral solution is that the 
former is given as an explicit and compact formula, which can be used to obtain the 



Power spectra Brownian motion with solid friction 



10 



propagator by inverse Laplace transform. The Laplace solution is also useful as it 
enables us to obtain the power spectrum of the system, in addition to the diffusion 
constant of the Brownian motion affected by solid friction. Both of these characteristics 
are easily accessible experimentally and so might be useful to compare the model 
with experiments involving noise and solid friction, such as those recently reported for 



example in 12 13 



Acknowledgment s 

H.T. is grateful for the support and hospitality of the National Institute of Theoretical 
Physics at the University of Stellenbosch, South Africa, where part of this work 
was written. The work of W.J. is partly supported by EPSRC through grant no. 
EP/H04812X/1. 



Appendix A. Laplace transform of the propagator 

We provide here some details of the derivation of the expressions |7|)-(|9]) for the Laplace 
transform P(x, s\x',0) of the propagator P(x, t\x', 0). 

The derivation starts with the Fokker-Planck equation ([5]), which for positive 
values for the initial condition x' > reads 

sP(x, s\x', 0) = ^{x-5- f)P(x, s\x', 0) + 

ax ax z 



sP(x, s\x', 0) = ^-(x + S- f)P(x, s\x', 0) + 

ax ax z 

d d 2 P 
sP(x,s\x',0) = —(x + S-f)P(x,s\x',0) + -—r, x'<x. (A.3) 
dx dx 1 

These equations must be solved with the following matching conditions at x = 0: 



x < 
< x < x' 

x' < X. 



(A.l) 
(A.2) 



(-S- f)P(x = 0-,s\x',0) + 



P(x = Q-,s\x',0) 
dP 



dx 



= p( x = 0+,s|x',0) 

= (S- f)P(x = 0+,s\x',0) + 



x=Q- 



dP 

dx 



x=0+ 



P{x = x' - 0, s\x', 0) = P(x = x' + 0, s\x\ 0) 



dP 

dx 



dP 

dx 



(A.4) 
(A.5) 
(A.6) 
(A.7) 



which results from the continuity of the probability current, as well as the usual 
decaying boundary conditions at x = ±oo. 

The equations ( A.l[ )-( A.3) have the form of the Hcrmite differential equation 

u" {z) + [zu(z))' + vu{z) = 0. (A.8) 

The solution can be written in terms of parabolic cylinder functions either as 
e~ 2 / 4 D„(z) or e~ z ^ 4 D l/ (—z). [+] In view of the asymptotic property of the parabolic 
cylinder function, 



D v {z)~3Te 



z — > oo, 



(A.9) 



+ These two solutions are linearly independent, i.e., they constitute a fundamental system if the 
index v is n ot an integer. For a complete account of parabolic cylinder functions, the reader may 
consult fl4l. 
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the solution of Eqs. (A.l) and (A.3) having vanishing currents at infinity can be written 
as 



P(x,s\x',0) = C-(x',s)e- {x - s - f)2/4 D^ s {~x + S + f) (A.10) 



for x < and 



P(x, s\x', 0) = C+(x', s) e -(*+ 5 -/) 2 /4 D _ s {x + S-f) (A.ll) 

for x > x'. On the other hand, the solution of Eq. (A. 2 1 is given by a linear combination 
of the two fundamental solutions: 

P(x, s\x', 0) = B+{x', s) e -(*+ 5 ~/) 2 /4 D^ s (x + 5 - /) 

+ B-(s,x , )e-( x+5 -V 2 / 4 D^ s (-x-6 + f), < x < x' . (A.12) 



The amplitudes B± and C± are determined by the matching conditions (A.4)-(A.7), 
which result in a set of inhomogeneous linear equations: 

C_(a;', s) e - {s+f)2/A D_ S (S + f) = B+(x', s) e - (s - f)2/A D_ s {5 - /) 

+ B_(s, x') e -( 5 -/) 2 / 4 D_ s (-6 + f) (A.13) 

-C.(x', s) e -^ 2 / 4 D- s -i((5 + /) - B+(x', s) e~^ 2 / 4 D_ s _i(<5 - /) 

- £_ (a/, s) e~^-rf/ 4 D_ s _ x (-«5 + /) (A.14) 
C+(x', s)£)_ s (x' + (5 - /) = B + (x', s)D_ s {x' + S-f) 

+ B-(x',s)D_ s (-x' -5 + f) (A.15) 
e W+B-f)*/A = s )^_ s+1 (x' + 5 - /) 

- s)D_ s+1 (x' + 6- f) 

+ B-(x',s)D- s+1 (-x' -S + f). (A.16) 

In writing these equations, we have used the following identities for the parabolic 
cylinder functions: 

ve^l A D v ^{z) = e- z2 '\zD v {z) - D v+1 {z)) 

= z e- 2 '/ 4 D v {z) + (e-* 2 / 4 D v (z))' (A.17) 

to evaluate and simplify the derivatives. 

Equations ( A.13)-(A.16| can be solved directly to find B± and C±. The difference 

(A.18) 



of Eqs. ( |AT5l ) and ( |A.16[ ) yields 

^le^'+^^D^x' + S-f) 



B_(x',s) 



20r 



if we take into account the product identity 

D v {z)D v -x{—z) + D v {-z)D v - X {z) = 



/2tt 



(A.19) 



which follows from the Wronskian of the fundamental system^] Then the difference of 
Eqs. \k.li\ and \k.\A\ yields 

C_(x',s)= g< (s,6,f)B_(x',s), (A.20) 



Equation (A.19i corrects a typo in Eq. (33a) of 14 . 
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<s _ ff/4 ( D_ s (-5 + f) D-s-ij-S + f) 

\ D- a {S- /) D- B _i(tf-/) 

_ (5+/)2/4 ( D- s (5 + f) D^jS + f) 



9<(s,5,f) = 

e -i.s+fY/A I 

\D_ s (5-f) ' D_ s _ x (<5-/) 

Summing Eqs. (A. 13 1 and ( A.14[ ), we also find 

B + (x',s) = g> (s,6J)B-(x',s) 

with 

£>_ s -i(-5 + /) £>_.(-* + /) 



9>(s,5J) 



D- fl -!(J + /) £>-,(<* + /) 
Ds(S-f) D-s-tf-f) 
D_ S (S + f) + + f) 
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(A.21) 



(A.22) 



(A.23) 



Finally, summing Eqs. (A. 15) and (A. 16), we find 



C + (x',s)=B + (x',s) + 



e (*'+«-/) 2 /4 #_ s („ x ' - ,5 + /). (A.24) 

V 27T 



The substitution of all these expressions back into Eqs. ( A.10)-(A.12) leads us to 
Eqs. 0"®- 

Appendix B. Asymptotic expansion 



The power spectrum, defined in Eq. ( 26 1 , is entirely determined by the ratio F s (5) of 
parabolic cylinder functions defined in Eq. ( 28 ) . Using the linear recurrence relation 
(A. 17), this ratio is found to obey the relation 

Fs{6) = 6 + (s + l)F s+1 (Sy 



(B.l) 



iteration of which leads to continued fraction expansions 11 



Asymptotic expressions for F S (S) for large values of S and s can be obtained from 
the contour integral representation of parabolic cylinder functions, 



D u (z) 



2/4 pc+ioo 



W2n 



-zt+S/2 t v dt 



(B.2) 



which is valid for all i,i/£C, and any integral contour satisfying | argi| < 7r/2 and 
c > 0. To this end, we introduce a large parameter N and change the variable of 
integration to u = t/\/N. This leads to 



F Ns (VNS) = 



c+zoo 



N(-Su+u 2 /2-slogu) 



da 



c— zoo 
c+zoo 



e N(-Su+u 2 /2-slogu) du 



Each of the integrals above has the Laplace form 

z Nf ^g(z)dz, 

which can be approximated following the saddlepoint method as 

2tt 



e Nf{z) g{z)dz ~ ±ie Nf ^g{z s ), 



Nf"(z s y 



(B.3) 



(B.4) 



(B.5) 
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where z s is the saddlepoint of the integral given by f'{z s ) = 0, and the sign depends 
on the orientation of the contour C, which is deformed to a steepest-descent contour 
passing through the saddlepoint. This approximation is valid to leading order in N 
provided g does not vanish at z s . 

In our case, we obtain the saddlepoints from 

<9/l \ s 

= — [ -du + -u 2 - slogu = -S + u , (B.6) 

du \ 2 J u 

which has two solutions: 

u x = l - (S + VS 2 + 4s) , u 2 = i (S - VS 2 + 4s) . (B.7) 

On the one hand, if 8 2 + 4s > 0, then u\ > > 1*2 and the steepest-descent contour 
passes through u\ only. In this case, much of the saddlepoint approximation cancels, 
leaving us with 

F Ns (y/NS) =— = -= ; B.8) 



On the other hand, if 6 2 + 4s < 0, then u\ and are complex conjugate to each other, 
and the steepest-descent contour passes through both saddlepoints. The crossover 
between these cases is given by coalescing saddlepoints; a uniform asymptotic expansion 



for this case can be obtained in terms of Airy functions 15 



A second case of interest is given when the argument of F is kept constant, i.e., 



c+ioo j 
c N(-8u+u 2 /2) aU 



F ^'^> ^r- (R9) 

00 u s 

In this case, there is only one saddlepoint at uq — 5, which collides with a singularity 
of the integrand when S approaches zero. For 5 > 0, we then find 

which matches the previous asymptotic evaluation. For d < 0, i.e., for Kramer's 
problem, the integral is dominated by the singularity at the origin, and its contribution 
gives 



/— -VN5 

F S (VNS) ~ (B.ll) 
s 



for s ^ 0, -1, -2, .... If s = -1, -2, . . ., we find instead F S (VNS) - 1/(VNS) 
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